FERMILAB-Pub-02/109-T 
hep-ph/0206058| 



Remark on the Theoretical Uncertainty in B°-B° Mixing 

Andreas S. Kronfeld 

Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia, Illinois 60510 

Sinead M. Ryan 

School of Mathematics, Trinity College, Dublin 2, Ireland 
June 5, 2002 

Abstract 

We re-examine the theoretical uncertainty in the Standard Model expression for 
mixing. We focus on lattice calculations of the ratio £, needed to relate the 
oscillation frequency of B®-B® mixing to the poorly known CKM element Vtd- We 
replace the usual linear chiral extrapolation with one that includes the logarithm 
that appears in chiral perturbation theory. We find a significant shift in the ratio £, 
from the conventional 1.15 ± 0.05 to £ = 1.32 ± 0.10. 

It is anticipated that the oscillation frequency of B®-B® mixing will be measured during 
Run 2 of the Tevatron [ |]]. It is thus timely to assess the measurement's impact on tests 
of the Cabibbo-Kobayashi-Maskawa (CKM) picture of flavor and CP violation. The CKM 
interpretation is limited by the poorly known hadronic matrix elements for B® «-> B® and 
SS <-> B® transitions. In this paper we re-examine lattice calculations of these matrix elements, 
focusing on the chiral extrapolation. We find that the range usually quoted is probably incorrect. 

In the Standard Model, the theoretical expression for the oscillation frequency is 

where q £ {d, s}, So is an Inami-Lim function, tjb is a short-distance QCD correction, and M. q 
is the hadronic matrix element for B® «-> B® transitions. In Eq. (jl]), the parentheses consists of 
accurately known quantities, and IV^V^ is the CKM factor. The hadronic matrix element 

M q = (B° q \[h^l - 7 5 MMl - 7 5 kK> (2) 



and r]B depend on the renormalization scheme, but the product risMq does not. The renor- 
malization-group invariant value of the short-distance factor is f/B = 0.55. 



One should keep in mind that non-Standard physics at short distances can modify Eq. ([!]). 
For convenience we shall couch the discussion as using Am q and the hadronic matrix element to 
determine | Vt q \ ■ The resulting value of | Vt q | can then be compared to other CKM determinations 
to test for deviations from the Standard Model. 

M. q must be computed with a non-perturbative method, such as lattice gauge theory. For 
historical reasons one usually writes 

M q = \mlj 2 B B Bq (3) 

and focuses on the decay constants f Bq and the bag parameters B Bq . But lattice QCD gives 
M. q directly (and f Bq separately from (0\try ^ q^®)) . The separation does, however, turn out 
to be useful, as we shall see below, when considering the dependence of f Bq and B Bq on the 
masses of the light quarks. 

At present the uncertainty in the matrix elements is large. A recent review [ || of lattice 
calculations quotes 

f Bs = 230 ± 30 MeV, B Bs = 1.34 ± 0.10, (4) 
f Bd = 198 ± 30 MeV, B Bd = 1.30 ± 0.12. (5) 

These estimates take into account the first (partially) unquenched calculations of f Bq [ |j| |||, [f| 
^, f?J, several quenched calculations of B Bq and preliminary results suggesting that B Bq changes 
little when the quenched approximation is removed [ g]. The raw Monte Carlo data in lattice 
calculations are generated with the light quark mass m q in the range 0.2-0.5 < m q /m s < 1, and 
the physical matrix elements are obtained by extrapolating m q to the down quark's mass m^. 
This method of reaching physically light quarks is called the chiral extrapolation, and it plays 
an important role in our analysis below. 

The frequency for B®-B^ mixing has been measured precisely, Arrid = 0.494 ± 0.007 ps _1 [ 
||]. With Eqs. ([!]) and (||) the uncertainty on \ Vtd\ is limited to 15% by fB d \jB Bd . The precision 
on \Vtd\ will not improve until better (unquenched) lattice calculations have been carried out. 
The frequency for B®-B® mixing is known to be high, Am s > 15 ps _1 [§. But details of the 
way Am s is extracted from the data mean that the first measurement will immediately have a 
precision at the percent level [ |J]. Thus, it is interesting to form the ratio 



Am s 
Am d 

where 



Vt 



ts 



v td 



2 ^ 2 , (6) 
m Bd 



f = (7) 

JB d B B d 

and use Eq. (^) to determine \Vtd\- The measurement uncertainties are (or soon will be) negli- 
gible. By CKM unitarity \Vt s \ = 1Kb I to good approximation. Thus, the error in \ V t d\ is 

s\v td \ = V / (M+W. (8) 

The uncertainty in |V^|, determined from semileptonic B decay, is also dominated by QCD, 
but it is only 2-4% and relatively well understood [ ^, 10, 11]. 



The conventional wisdom, coming from several reviews of lattice B physics, is that <5£ is 
small. Based on such endorsement, recent efforts to fit a wide range of precisely measured 
flavor observables have used £ = 1.14 ± 0.03 ± 0.05 [ || or £ = 1.16 ± 0.03 ± 0.05 [ The 



2 



second error bar is meant to reflect the uncertainty from the quenched approximation; the first 
covers all other sources of uncertainty in lattice calculations. Central values in this range are 
reproduced by many quenched, and some unquenched, calculations. 

Such a small error is, however, not universally accepted in the lattice community. Booth [ 
|l4|| , noting that chiral logarithms in the quenched approximation differ strikingly from those of 
QCD, predicted that £ in QCD would be 0.15-0.28 larger than in the quenched approximation. 
Sharpe and Zhang [ [T!|, with a similar point of view, reckoned that <5(£ — l)/(£ — 1) could 
be 100%. Bernard, Blum and Soni [|l6| studied two different ways of carrying out the analysis. 
Treating /s s //s d and BbJBb a separately (as usual), they found £ = 1.17 db 0.02 j^gg; treating 
instead M s /M d directly, they found £ = 1.30 ± 0.04±g;^. (In Ref. [|l| the second error comes 
from studying the lattice spacing dependence; the difference is significant source of concern [ 
p7|.) Finally, the JLQCD collaboration studied the effect of the chiral log in lattice calculations 
with two light flavors, finding that the extrapolated value of £ could change significantly [|(|. 

At first glance, could well be smaller than <J/b,//b 9 - £ is a ratio of similar quantities, 
so, in numerical lattice calculations, most of the Monte Carlo statistical fluctuations do cancel. 
Similarly, the short-distance normalization factor of the lattice operator also cancels. But one is 
still left with a multi-scale problem, with the heavy quark mass rrij,, the QCD scale Aqcd and the 
range of light quark masses from m s down to m^. Because the numerator and denominator of £ 
are the same, except for the light quark, one may expect £ to be insensitive to the heavy-quark 
and QCD scales, but not to scales between m s and md- 

Let us examine the uncertainties associated with each scale in more detail. Heavy-quark 
corrections to £ are suppressed by (m s — m^/rrib ~ 2%. In lattice calculations, one should 
also worry about discretization effects of the heavy quark, because m^a ~ 1. There are several 



ways to handle this problem and some debate over the best method [IS]. But the various 
discretizations yield consistent results for fB 3 /fB d and Bb 3 IBb a - Thus, we conclude that 
errors from the short distance scales are under control. 

Next let us consider Aqcd- Implicit in the quenched approximation (also called the valence 
approximation) is that the omitted sea quarks are compensated by a shift in the bare gauge 
coupling [ P^9[ . This treats light-quark vacuum polarization in a dielectric approximation. Such 
approximations can be accurate when looking at a narrow range of scales. In the case at hand, 
that means that ratios of decay constants or bag parameters could be accurate as long as all 
quark masses are not too different. Thus, it is plausible that the quenched approximation accu- 
rately determines the slope of £, viewed as a function of r = m q /m s , when r ~ 1. Unquenched 
calculations [ ||, H g, |||, ]7j do not contradict this expectation. These calculations, and the 
justification of the quenched approximation [ 19 1, suggest that the scale Aqcd is also under 
control. 

That leaves us with contributions to £ from the long distances between l/m s and 1/uid- 
Here the quenched approximation is known to break down [ [ll| , and it is not obvious that 
the quenching error could be as small as 5%. One must take a careful look at how the chiral 
extrapolation is done, and consider what methods of extrapolation are reliable. 

The correct framework to discuss the long-distance behavior of QCD, and the chiral extrap- 
olation in particular, is chiral perturbation theory. We neglect 1/m corrections and write 

^/m-j Bq = $ [i + Af q ] , (9) 

B Bq =B[1 + AB q ] , (10) 

where $ and B are independent of both heavy and light quark masses, and Af q and AB q 
denote the (one-loop) contribution of the light meson cloud. The "chiral logarithms" reside in 



3 



Af q and AB q . 



Neglecting isospin breaking, the one-loop corrections to the decay constants are [ ^0], ^lj, [L4|, 



m 



Af s 



Afd 



(4tt/) 2 



l + 3g 2 
(4vr/) 2 



[ml- In {m 2 K /ft 2 ) + \m 2 In (m 2 //i 2 )] 



+ (m| + imj) ftifi) + (m| - ±m 2 ) / 2 (/x), (11) 



[§m 2 In (m 2 /^ 2 ) + |m|f In [m 2 K /n 2 ) + ^m 2 In (m 2 /^ 2 )] 
+ + \ml) fi(fi) + im 2 / 2 (M), 



(12) 



and to the bag parameters 
1 



AS, 



' 3g 2 2 2 i / 2/2 
(4^)2- 3 m v ln K/^ 



+ 



/ 2 , 1 2^ 

(m x + 



+ (m 2 K 



>m„ 



A5,, 



l-3g 2 

(47T/) 2 



[±m 2 ln(m 2 //z 2 ) + ±m 2 In (m 2 //i 2 )] 



S 2 0*), (13) 



(14) 



where / and g are (the chiral limit of) the light pseudoscalar decay constant and B-B*-ir 
coupling. The "low-energy" constants /i(/u) and Bi(fi) encode QCD dynamics from distances 
shorter than fjf 1 , whereas the logarithms are long-distance properties of QCD, constrained by 
chiral symmetry. The dependence on fx cancels in the total. 

It is convenient to look separately at the fs and \/Bb factors in £. The chiral logarithm 
in the y/Bs factor could be small because it is multiplied by 1 — 3g 2 . On the other hand, 
the chiral logarithm in the /b factor could be significant, because it is multiplied by 1 + 3g 2 . 
Consequently, we focus on 

tf = f B JfB d (15) 

and study its chiral extrapolation. Our strategy is to use lattice calculations as an (indirect) way 
of determining the low-energy constants, and then we reconstitute £/. Repeating our analysis 
for the chiral extrapolation of £b = \J Bb s I B^ d verifies that £b has a small effect. 
Combining Eqs. (11) and (|j~2|), the first non-trivial order in the chiral expansion is 



e/-l = (m|-m 2 )/ 2 (/x) 



l + 3g 2 

(47T/) 2 



[\m 2 K \n{m 2 K /n 2 ) + \m\\v,{m\l \i 2 



§m 2 ln(m 2 //i 2 )] . 

(16) 

All lattice estimates of £ are obtained not at physical light meson masses, but by chiral extrap- 
olation. Therefore, we use Gell-Mann-Okubo formulae to replace the meson masses with 



in 



qqi 
2 



(m 2 s + m 2 ? )/2, 
(2m 2 , + m 2 )/3. 



(17) 
(18) 
(19) 



Varying the light quark mass changes m 2 oc m q . Lattice calculations typically take m 2 not 



too different from m 2 , so we write m 



rm: 



Then 



£/(r)-l = m 2 s (l 



l + 3g 2 

(47T/) 2 



12 



ln(m 2 s /^ 2 ) + l(r) 



(20) 
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where 



l(r) 



1 



1 + r 



In 



1 + r 



+ 



12 



In 



2 + r 



3r 



ln(r) 



(21) 



The function x(r) = (1 — r)l{r) contains the chiral logarithms. It is plotted in Fig. |l[ The 
curvature over 0.5 < r < 1.0 is too small to be resolved when there are percent-level statistical 
uncertainties on £j. But once r <C 1, which is appropriate for the down quark with r^ ~ 
1/25, the curvature required by the chiral log has a significant effect. Fig. || shows this effect, 
comparing the conventional linear chiral extrapolation with Eq. (p0|), for f^i^) in the range 
coming from Eq. (|2~3|), below. 

When £ is calculated in lattice gauge theory, the range of r is restricted to r < 1 but r <jL 1. 
Usually, it is fit to a straight line 



t f (r)-\ = (l-r)S f 



(22) 



and similarly £fj(r) — 1 = (1 — r)SB- Usually one assumes this linear extrapolation holds down 
to the chiral limit, quoting £ = [1 + (1 — 7-^)5/] [1 + |(1 — r^Ss]- The chiral log says, however, 
that this procedure is not trustworthy. It has been employed because there was, until recently, 
no independent reliable estimate of the B-B*-ir coupling g 2 in the coefficient of the chiral log. 

The CLEO collaboration has recently measured the width of the D* meson, which yields a 
value for the D-D*-ir coupling [ ^J. Heavy-quark symmetry suggests that the B-B*-ir coupling 
is nearly the same. On this basis, we shall set g 2 = 0.35, although below we allow for 20% 



deviations. With g 2 = 0.35, the chiral log in is truly small, because 1 - 3g 2 = -0.05, but 
the chiral log in £j is multiplied with 1 + 3g 2 = +2.05. 

With this handle on g 2 , we can interpret the lattice results for Sj as a calculation of /2 (/•*)■ 
We assume the linear fit given by Eq. ( [22] ) makes sense around r = r$ ~ 1, even though we do 
not trust it when r«l, So, at ro we set the right-hand side of Eq. ( |20| ) equal to the right-hand 
side of Eq. (H) and find 



mtskhiv) =Sf + m 



Then, inserting this result into Eq. (|2 



2 l + 3g 2 
(4vr/)2 



— ln(m 2 Jn 2 ) +/(r ) 



(23) 



fr(r) - 1 = (1 - r) Is, + ™L^|C [J(n>) - «r)]\ ■ 

To evaluate the right-hand side, one needs estimates of /, g 2 and Sf. We use / 
and g 2 = 0.35. In addition, we take [|2| 

(1 -r d )S f = 0.15 ±0.05 



(24) 
130 MeV 

(25) 



which brackets many quenched calculations (for which there is a lot of experience and repro- 
ducibility) as well as less well-developed unquenched calculations .[] 

Once we have made the Ansatz to use the slope from lattice QCD to determine the low- 
energy constant via Eq. fl23|), another source of uncertainty is the choice of ro. Fig. || shows 
the result from Eq. (|24|) for the physical value Cf{ r d)^ as a function of ro from to 1.5. (The 
lower end is not natural, but recovers the conventional result; the upper end 1.5 is where this 
order of chiral perturbation theory is less trustworthy.) Since the typical range of fits leading 



1 In fact, some "unquenched" calculations only have nf = 2. 
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Figure 1: Plot of the chiral logarithm x( r ) as the mass ratio r = rn^ q /m 2 ss = m q /m s is varied, 
compared with a straight line fit for 0.5 < r < 1.0. The difference between the curve and the 
fit is shown in the inset. 
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to Eq. (|2q) is 0.5 < r < 1.0, we choose tq in this range and use Fig. || to obtain 

i f = 1.32 ±0.08. (26) 

With separation scale /j, = 1 GeV, £j — 1 receives nearly equal contributions from the low-energy 
constant (0.159) and the chiral log (0.165).Q 

We have carried out a similar analysis for £g and also allowed for a ±20% range on g 2 . 
(See the appendix for details.) The chiral logs in £y and is pull in opposite directions, so the 
resulting £ = £/£s is insensitive to g 2 : 

£ = 1.32 ±0.10, (27) 
which is quite different from the range usually used in CKM fits, although it agrees with 



qualitative discussions of chiral logs [ |bj, |15[| , a direct analysis of M. s /M.d [ |16|] , and chiral log 
fits to preliminary unquenched calculations [ ||]. The shift in central value from 1.15 to 1.32 can 
be thought of as a correction to the quenched approximation: mature unquenched calculations 
will certainly see the curvature required by the chiral log. 

Because our result is so different than the conventional one, let us stress the differences in 
methodology. Usually £ is obtained via a linear chiral extrapolation, although chiral log fits 
have been tried in Ref. [ |(|. We have relied completely on the functional form predicted by 
chiral perturbation theory. It is difficult to determine the coefficient of the chiral logs directly 
from the lattice calculation. We have circumvented this obstacle by using the D* width [ |22j , 
which, with heavy-quark symmetry, implies g 2 = 0.35. The uncertainty in Eq. (|27| ) is larger 

2 Loops with excited B** mesons are expected to contribute significantly to £j> [ ^i), but the ensuing r 
dependence is well described by linear extrapolation, so it is accurate to lump them into (1 — r)/a(/i). 
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than in many other papers, mostly because we have assigned ±0.05 instead of ±0.03 uncertainty 
to the lattice calculations. On the other hand, we have also not done a complete error analysis: 
for example, we have neglected uncertainties from higher orders in the chiral expansion. 

One could easily reduce the theoretical uncertainty in B°-B° mixing by carrying out lattice 
calculations designed to determine the low-energy constants in Eqs. (|IT|)~(|l4|). If one takes 
closely-spaced values of the light quark mass, even if close to the strange mass, one can compute 
the derivative d^/dr. If one is willing to take g 2 from experiment, these derivatives give ^(/-t) 
and £?2(m); an d one can proceed to determine £ for physically light quark masses. The same 
procedure could be applied to /r and Bs q although now one must also compute and 
Bi(n), and also cope with further low-energy constants in the l/m& corrections [ ^]. Chiral 
extrapolations with chiral logs may well change /g from the estimates in Eqs. (||) and (|5|) in 
the same way they changed £/. 

From a (lattice) purist's point of view it may be unsatisfactory to take g 2 from experiment. 
In the long run it will, however, be possible to solidify our knowledge of g 2 (in the B system) 
through lattice calculations and other applications of chiral perturbation theory to B physics. 
To relate the very precise measurements to the CKM matrix, the combination of phenomenology 
for g 2 and lattice calculation for the low-energy constant is very satisfactory, especially since we 
find that £ varies by less than 2% when g 2 is varied by 20%. Fig. || shows how the combination 
of sin 2/3 and Am s / Arrid work together to constrain the apex of the unitarity triangle. We 
take sin 2/3 = 0.79 ± 0.10 from averaging CDF [ 27|, BaBar [ 28] and Belle [|2j| measurements. 
For illustration we take Am s = 20 ps _1 , and compare £ = 1.15 ± 0.05 (conventional wisdom) 
with £ = 1.32 ± 0.05 [Eq. (|27| ) with error halved]. With a larger value of £ the mixing side is 



longer, scaling like £/y/Am s . By the same token, our result suggests that the Standard-Model 
prediction for Am s (16-19 ps _1 [|l2|]) should be increased, perhaps by 25-35%. 
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Figure 4: Sketch of constraints on the apex of the unitarity triangle with sin 2(3 = 0.79 ± 0.10, 
Am s = 20 ps" 1 and £ = 1.32 ± 0.05 or 1.15 ± 0.05. 
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Appendix: Analysis including £g and varying g 2 

Let £ B = BB s /BB d , with linear chiral extrapolation 

C|(r)-l = (l-r)5 B . 
Then, eliminating ^(/x) in the same way as /2(m) m £/> 



el(r) 



;i-r)<j5 B + mL^|^ 



where 



k(0 



1 



1 



2 + r 



In 



2 + r 



[^(ro)-^(r)] 



ln(r) 



To evaluate the right-hand side, we take [ || 

S b = 0.00 ±0.05. 



(28) 

(29) 
(30) 

(31) 



Then we find £ B = 0.998 ± 0.025. 

In the main analysis, we have used g 2 = 0.35, which assumes that the B-B*-n and D-D*-ir 
are the same. Repeating the analysis with g 2 = 0.20 and 0.50 we find the results in Table [l|. 
Although the chiral extrapolation of £g is no longer completely insignificant, and £y changes a 
little, the result for £ is very stable. 



Table 1: Comparison of chiral extrapolations for £f, £g and £ for three values of the B-B*-tt 
coupling g 2 = 0.20, 0.35, 0.50. 
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£s 


e 


0.20 


1.29±0.08 


1.01±0.03 


1.30±0.09 


0.35 


1.32±0.08 


1.00±0.02 


1.32±0.09 


0.50 


1.36±0.09 


0.98±0.02 


1.34±0.09 
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